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In this paper we derive the dispersion relation of the surface waves at the interfaces between Mott 
insulating and Superfluid phases for confined Bose-Einstein Condensates. We then calculate their 
contribution to the heat capacity of the system and show how its low temperature scaling allows an 
experimental test of the existence and properties of Mott Insulator-Superfluid domains. 
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Following the experimental observation of Bose- 
Einstein Condensation (BEC) pj a lot of attention has 
been given to the investigation of the possible quantum 
phase transitions such a system can undergo. Within the 
framework of the Bose-Hubbard model, Fisher et al. 
proposed the zero temperature phase diagram of an in- 
finite uniform system on a lattice as a function of the 
typical hopping kinetic energy J, the atomic on-site re- 
pulsive interaction U and their chemical potential [i. In 
the (J/U,fi/U) plane they found lobe-like phase bound- 
aries between Mott-Insulating (MI) and Superfluid (SF) 
phases for small and large values of J/U, respectively. 
Recent experiments on BEC in a periodic optical lat- 
tice showed the expected MI-SF transition 0, 0] driven 
by variations of the system parameters, while theoretical 
investigations focused on the restoration of macroscopic 
phase coherence after abrupt parametric changes 0, |g . 
Real experiments are however always performed on fi- 
nite systems subject to an external confinement, yield- 
ing an effective space-dependent chemical potential. The 
consequent modulation in the local boson density can 
produce spatial-domains characterized by different equi- 
librium phases (MI or SF). Indeed, it has been recently 
proposed that the observed MI-SF crossover is to be in- 
terpreted as a relative widening/shrinking of the different 
phase domains rather than as a pure infinite long-range 
transition [tJ. 

Recent numerical calculations investigated the equilib- 
rium state of the confined BEC and confirmed the exis- 
tence of the above-mentioned MI-SF domains 0, El • 
However, their experimental signatures are still an open 
issue. Moreover, it is clear that the interfaces between the 
different phases can play an important role in the ther- 
modynamic and time-evolution properties of the system 
under external perturbations. 

In this paper we consider the interfaces between MI and 
SF phases and evaluate the dispersion of their surface 
wave excitations. We show that the low temperature 
heat capacity of the system is crucially affected by the 
existence and properties of MI-SF domains, allowing for 



an experimental way to detect them. 

A BEC at zero temperature in a periodic optical lattice 
can be described by the Bose-Hubbard Hamiltonian 



H 



J Y, ( a k- + h - c ) + J2 + V ^> a K + 



u 



(1) 



with a i the annihilation operator for a boson at site i 
and (ij) indicates nearest-neighbor hopping. The single- 
particle kinetic term J can be tuned by the intensity of 
the light generating the optical lattice, the on-site repul- 
sion U by the use of Feshbach resonances while the chem- 
ical potential /i depends on the particle number and V i 
describes an external confining potential. 
For the uniform case V i = 0, the system described by 
the Hamiltonian Q undergoes quantum phase transi- 
tions between MI phases and SF ones, driven by the two 
parameters J/U and fi/U. In the Mi-phases the number 
of particles per-site is fixed to a positive integer, produc- 
ing a constant density, in the continuum limit. 
Under the influence of a confinement, the effective chem- 
ical potential [i — is position dependent and the sam- 
ple can show MI phases and SF ones in different do- 
mains. We will consider an effective two-dimensional sys- 
tem in the (x, z) plane with a one-dimensional confine- 
ment along the z-direction and translational invariance 
along the x-axis. A realistic cylindrically symmetric con- 
finement reduces to our model by remapping the radial 
direction into the z-axis and considering the x direction 
as closed per periodicity. In principle many successive 
alternated MI and SF phases can be formed, according 
to the particle density, but we will mainly concentrate 
on one SF domain "guarded" by two neighboring MI 
ones. The confinement is such to produce a Mi-phase 
with fixed integer site occupation n M for z < and the 
next MI with occupation n M + 1 for z > L, with a SF- 
phase with variable density for < z < L. If we use, 
for simplicity, a square lattice with typical lattice size I, 
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in the continuum the two MI have densities p M = n M /l 2 
and p M + Ap = (n M + l)// 2 . The equilibrium SF density 
p{z) varies continuously between these values, 



p(z) = p M + f(z) Ap 



(2) 



where f(z) is determined in principle by solving the 
Gross-Pitayevski equation in presence of the external 
confinement and is such that /(0) = and f(L) = 1. 

We imagine to impose an infinitesimal displacement 
£(x, t) oc exp[i(fca; — cot)], £ <c L, along the z-direction to 
the interface at z — while that at z — L does not move. 
Our task will be to determine the change in the "potential 
energy" associated to this perturbation via the surface 
tension of the interface, as well as the "kinetic energy" 
due to the induced SF velocity. Up to quadratic terms 
in £ and £, for any k, the former contribution will have a 
generic form 1/2 Mu; 2 (k)( 2 and the latter 1/2 MC, 2 with 
M a generalized mass, whence the surface waves disper- 
sion Lo(k). 

This problem is similar to the determination of the dis- 
persion for crystallization waves at a quantum interface 
between solid and liquid He 4 |llj . However, the continu- 
ity of the density profile and the presence of the cutoff 
length L bring important new features into our case. 

Let us first calculate the induced change in the poten- 
tial energy of the interface 



U 



dl a(n) 



(3) 



where dl is the unit length element along the surface and 
a(n) is the surface energy per unit length as a function of 
the normal n to the interface. In the k = case, the rigid 
shift of the interface away from its equilibrium location 
will correspond to an energy cost 1/2 A £| 2 , whose har- 
monic confinement frequency A will be analyzed later. 
For finite-fc and small £ we can expand n z = 1 and 



n x = d x (, so that dl = dx 1 + 1/2 (d x ()' 
face tension term undergoes the expansion 



The sur- 



a(n) = a(z) + 



da 
dn^ 



d x C + 



1 d 2 a 

2 dri 2 



(d x Cf (4) 



so that the variation per unit-length in the potential en- 
ergy of the interface due to the modulation is 



AU 



A 



a(z) 



di 



k 



icr 



(•5) 



where we used the relation d x ( = ikQ. 

Now let us concentrate on the kinetic energy associ- 
ated to the superfluid velocity v induced by the surface 
displacement. In order to determine the velocity profile 
we consider the continuity equation 

d t P\ oc {x,z,t) + y -[p loc {x,z,t)w{x,z,t)\ =0 (6) 



with p loc (x,z,t) the local density. In the k — case, the 
density as well as the velocity do not depend on x and 
the only solution for v (the z-component of v) yielding a 
vanishing equilibrium current at z > L is 



v(z,t) 



1 



dyd t p loc (y,t) 



(J) 



The consequent uniform velocity in the MI phase at z < 
is needed to compensate for the loss (or addition) of 
particles implied by the variation in the density profile. 
The related current running in the insulating domain is 
therefore associated to a very slow relaxation process. 
In the finite-fc regime the two components of the velocity 
v = (v x ,v z ), show up. The SF constraint V x v(x, z, t) = 
is accounted for by defining v as a potential gradient, 
v(x, z, t) = Vip(x, z, t). 

In the limit of small oscillations, we write the density 
in terms of its equilibrium profile and a small induced 
fluctuation, p loc (a;,z,i) = p(z) + p{x,z,i). Thus, p and 
the related velocity v are the small quantities around 
which to expand the continuity equation JSJ to obtain 

^-*V + a,[ios^ + | = o . (8) 

The generic features of the surface waves we consider 
can be grasped by a model where we assume a linear de- 
pendence of the equilibrium density on the z coordinate, 
i.e. in @ we substitute f(z) = z/L. Such a description 
is relevant for either a linear external confinement or for 
the regions where a generic confinement can be effectively 
linearized. In any case, corrections to this model due to 
the curvature of the realistic confinement should not al- 
ter the main qualitative effects we now investigate. 
Eq. © can be solved in the limit pin/Ap ^> 1, (i.e. 
large site-occupation n M ^> 1). Then, i9 z [logp] = 2q, 
having defined q = Ap/ (2Lpyi), with the dimensions of 
a wave- number, and gL < 1. In the small occupation 
limit (always at least n M > 1) this approximation is in- 
validated only close to z — L and no major qualitative 
features will change. 

In order to specify the inhomogeneous term p/p, we con- 
sider the dynamics via the linearized Newton's law 

mp(z) 9 t v(x, z,t) H — Vp(x, z, t) = (9) 

K 

with m the mass of the single bosonic atom and n — 
dp/ dp the compressibility of the superfluid, assumed con- 
stant for simplicity. Eq. , together with the linearized 
continuity equation, lead to the familiar wave equation 



d 2 p{x, z, t) - c 2 V 2 p(x, z, t) = 



(10) 



yielding bulk sound waves with group velocity c = 
(mc)" 1 / 2 = \Zti Pm 1 2 U /m^ 12 [l^. The solution to the 
wave equation is of the form p oc exp[i(/ca; — cut) ± pz] 



3 



with p 2 = fc 2 (l — ui 2 /c 2 k 2 ), while the interface modula- 
tion corresponds to a boundary condition p(x, z = 0, t) = 
—2qp M ((x, t) and we still preserve p(x, z — L,t) = 
The solution to eq. i|l(J|) yields 

sinh[p(z — L)\ 



P -^2 QC ... 
p s\xm[pL\ 

with the familiar large-p exponential decay 



S -2qe~ pz C 



as well as the low-p expansion 



(ii) 



(12) 



(13) 



determined by the cutoff length L being smaller than the 
natural 1/p decay length. Alternatively, Eq. 113fl can 
be obtained using the adiabatic approximation for slow 
low-fc modes, replacing f{z) in (J2J) with f ad (x,z,t) — 

(*-C(ar, *))/(£ -C(M)). 

In analogy, in the theory of the crystallization waves for 
He the induced velocity and density modulation in the 
infinite SF region decay as exp[— kz] and are proportional 
to the difference in the solid and liquid uniform equilib- 
rium densities. In our case the density difference is rather 
replaced by the "slope" q since, if q = (i.e. Ap — 0) no 
fluid would move due to the displacement of the interface. 
Equation (JSJ thus has the general solution 



i/>(x, z, t) - 
(p 2 



c <y( k ) ex P [-QA k ) z \ 



(14) 



=± 



k 2 ) sinh[p(£ — z)] + 2qp cosh[p(L — z)} 



sinh[pi] ((p 2 — k 2 ) 2 

,2/„2\V2 



Ap 2 q 2 ) 

with Q a (k) = q l + cr(l + fc 2 /g 2 ) ±/i . In |H|) the co- 
efficients c±(k) are determined by imposing a vanishing 
v z = d z ip at both z — and z = L. In the finitc- 
k limit this corresponds to oscillating superfiuid motion 
with non- vanishing v x (z = 0) and v x (z = L). As a re- 
sult, if we concentrate on slow surface waves (p ~ k), the 
condition qL <C 1 yields the expansions 

1 i 



lfcL>l 



(qz 1 



-»C- f-(l + fez) e 



— fcz 



* lfcL» 



— fc2 



with the factor i meaning a 7r/2 shift between the in- 
terface modulation £ and the velocity v x . The velocity 
is proportional to ( and becomes very large for small fc. 
Thus, £ has to be thought of as infinitesimally small to 
prevent the SF from exceeding the critical velocity after 
which superfluidity is destroyed. 

The kinetic energy per unit length of the interface can 
therefore be calculated as 



K 



1 



dz p{z) 



G(fcL) 



(15) 



with G(x) a smooth function such that G(x «C 1) — 1 
and G(x 3> 1) ~ 3/(2x). Its universal properties are 
captured by the choice G{x) = [1 + (2x/3) 2 }- 1 ^ 2 . 

From the potential energy (J5J and the kinetic term 115|) 
we derive the dispersion of the surface waves 



Uj(k) ~ (3k 




1/2 



I 1/2 



with /3 = [A/ (p M mq 
low-fc linear regime and fc A = 



(16) 

the group velocity in the 

r / »2 \-|l/2 

_ A / («(*) + • 

For k > 1/L the interface at z = i plays no major role 
and we recover the fc 3 / 2 dispersion common to crystal- 
lization and capillary waves. 

We now proceed to calculate the thermodynamic prop- 
erties of the system, starting from the energy due to the 
occupation of the surface waves at temperature T, 



E = Ti 2 



/ dwtV 1D) (M- 

Jo > 



1 



,Ru/K B T _ l 



(17) 



with K B the Boltzmann constant and p( 1D ^ (hu>) the one- 
dimensional density of states. In the regime where the 
modes have a dispersion u> ~ jk a this is 



.(ID) 



(hco) 



2ir 



a (Jvy) 



l/c 



l/a-l 



(18) 



with L x the interface length (L x 3> L). The resulting 
contribution to the heat capacity Cv 10 ' = d T E is 



(1D) _ r(l + l/a)Z(l + l/q) 

— i / . 



2tt 



l/a 



1/'- 



(19) 

with r and Z the Euler and Riemann functions. This 
has to be compared with the bulk contribution due to 
the SF modes with ui = c (fc 2 + fc 2 ) 1//2 . The wavenumber 
quantization along z gives fc z = rnr/L (n 6 N), yielding 
several ID branches with energy split E g = hcir/L for 
vanishing k x . For K B T ^> E 1 the bulk behaves as two- 
dimensional with 



C (2D) 



L X L 6Z(3) 



2tt 



h 2 c 2 



KlT 2 



(20) 



For K B T < E only the n = branch is involved in the 
heat capacity, yielding a contribution analogous to l|19|) 
with a = 1 and 7 replaced by c. 

Thus, in the low-temperature regime, the heat capacity 
is dominated by the ID contribution and is linear in T. 

Before analyzing the implications of the results above, 
we estimate the typical parameters in our calculations, 
starting with L. The phase diagram for the infinite sys- 
tem shows that the transition between a MI with on- 
site occupation n and the adjacent one with n + 1, for 



4 



a given ratio J/U, implies a maximum jump fi/U = 1, 
independent on n. Thus, the effective chemical poten- 
tial n(z) = /t — V(z) undergoes a jump of the order of 
U between neighboring Mott-phases. If we describe the 
external confinement with the linear function V(z) = r\z 
this implies that the typical separation between one MI 
and the next is L ~ U /rj. 

As far as A is concerned, we estimate the variation of 
energy due to a small static density fluctuation p(r) 



6U =\J drdr'p(r)[/(r-r>(r') 



(21) 



with the on-site boson interaction U(r—r') ~ Ul 2 S(r— r'). 
In the k = case, via the adiabatic density fluctuation 
used in and from eq. © we get A = U/(3l 2 L). 
Last, we are left with a{z) while the contribution 
d 2 a/dn 2 is neglected as for quantum rough surfaces. We 
consider a straight interface as an ideal line separating 
lattice sites in the MI phase from those in the SF one. 
Neighboring sites in opposite phases are called " links" , 
and the surface energy per unit length is the product 
of the density of links (1/2) times the energy per link J 
given by the hopping term in leading to a(z) = J /I. 
Out of our previous definitions and estimates we finally 
get k A = [rj/MJ] 1 ^ 2 as well as f3 — ^2/3-kc. 
In the limit of large occupations n M ^> 1 the critical 
ratio (J/27) crit needed to drive the system into the SF- 
phase scales as l/n M . Thus, in particular, we also need 
U > n M J in order to have interfaces at all. In parallel, if 
the potential energy difference between neighboring sites 
is equal to the MI gap U, resonant tunnelling destroys 
the insulating phase Thus, in order to preserve the 
Mott phases, the confinement potential has to fulfill the 
condition r\l < U. These inequalities produce k A L 3> 1. 
In the limit J/U <C 1 and if the potential is steep, it is 
possible to achieve the regime L < I. Here the SF rings 
disappear, yielding neighboring MI domains whose sta- 
bility was analyzed recently 13]. In this case the ground 
state is everywhere incompressible with energy gap U and 
an activated heat capacity. 

As a result, while the bulk heat capacity is linear in T 
up to E g /K B and grows as T 2 for high temperatures, the 
contribution due to the ID surface waves is linear in T 
up to h(3/LK B , it scales as T 2 / 3 between h(3/LK B and 
Tiu{k A )/K B to finally grow as T 2 / 5 for K B T > huj{k A ). 
Thus, a direct measurement of the low-temperature heat 
capacity of the system will probe both the existence and 
the typical dimensions of the SF rings. Indeed, the very 
fact that SF domains exist introduces the cutoff length L 
responsible for the transverse discretization of the bulk 
spectrum and the low-fc linear dispersion of the surface 
waves. As a result, a linear- T dependence of C v signals 
the existence of MI-SF domains. 

Otherwise, if the atomic cloud is in a pure SF phase we 
expect a bulk T 2 scaling down to very low temperatures. 
Even in this case there is an extremely small quantization 



energy E associated to the lowest harmonic mode with 
half wavelength equal to the total cloud diameter. For 
K B T < E the system behaves as zero-dimensional with 
an exponentially activated heat capacity. 
Further informations about the number and typical size 
of the SF domains could be extracted from the slope of 
the linear-T part of C v . Indeed, the total heat capacity 
is given by summing the bulk and interface terms up to 
the number of SF domains N, 



cl lin) 



T(2) Z(2) K B T 
2~7rft 



N 

E 

71=1 



1 



fa 



(22) 

Since the sound velocities f3 n and c n in the n-th domain 
can be controlled via the external parameters, the slope 
of Cy yields informations about N and the typical 
length of the SF rings. Using typical experimental pa- 
rameters indicated in |4| we can estimate the temperature 
T g = E g /K Bl below which we expect the ID behaviour, 
to be of the order of T g ~ 200 nK, offering a wide window 
for experimental detection of this effect. 

In conclusion, after calculating the dispersion of the 
collective excitations localized at the MI-SF interfaces, 
we deduced their contribution to the thermodynamic 
properties of a confined 2D Bose system on a lattice. 
A direct measurement of a linear low-T heat capacity 
would signal the existence of MI-SF domains and its 
slope would tell about their typical size via the external 
control on J, U and the confining potential. At higher 
temperatures, in presence of domains, a crossover to 
other power- laws (C v - T 2 / 3 , T 2 / 5 ) is to be expected, 
before the SF bulk T 2 scaling dominates the thermo- 
dynamic properties of the system. The latter is the 
only leading contribution to C v in the absence of MI-SF 
domains, for a pure superfluid phase. 
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